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1.  Introduction: 


Let  B  be  a  p*p  symmetric  matrix  having  the  Wishart  distribution 
(1.1)  W  (B 1 1 1 f ) dB  =  C  f | B |  (f"P-1)/2  e“1/2  trB  dB  , 


where 

d.2)  c'i  -  2^  n 

pt  i-i  '  2 

and  dB  stands  for  the  product  of  the  differentials  of  the  p(p+l)/2 
distinct  elements  of  B.  Let  x  and  £  be  two  vector  variables  of  p 
components,  distributed  independently  of  B,  and  also  independently  of 
each  other,  as 


(1.3) 

and 

(1.4) 


1  -1/2  x'x  , 

e  - dx  , 


(277) 


P/2 


(2tt) 


A-  e'1/2  dZ 


P/2 


respectively.  While  considering  the  problem  of  multivariate  statistical 
outliers,  Wilks  (1963)  used  statistics  of  the  type, 
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(1.5)  r  ■  |b  +  |/|B  +  xx'  +  |  ,  s  ■  |B  +  xx'  | / j B  +  jot'  +  |  . 

He  has  remarked  that  the  exact  distribution  (joint)  of  r  and  s  is  compli¬ 
cated  and  has  given  the  expected  values,  variances  and  covariance  of  r 
and  s.  Unfortunately,  his  expressions  for  the  variance  and  covariance 
are  in  error.  The  purpose  of  this  note  is  to  derive  the  exact  joint 
distribution  of  r  and  s  and  to  give  correct  expressions  for  the  moments. 


(2.5) 


1  -  v'v 


|  B+xx'+yy 1 


Also  observe  that  in  (2.3) 


(2.6) 


I  -  uu'  -  vv_'  |  ■  (1  -  u'u)  (1  -  v^’v)  -  (u'v.) 


=  rs  -  (u'v) 


In  (2.3)  ,  transform  from  v  to  w  =  [w^  ,  w^  ,  . ..,  w  ] '  ,  by  an  orthogonal 


transformation 


(2.7)  w  =  Lv  , 


where 


L  is  a  p*p  orthogonal  matrix,  whose  last  row  is  u'//u 'u  .  The 
Jacobian  of  this  transformation  is  I L 1  *=  1  and  v'v  =  w'w  =  1-s  .  Also 


(2.8)  u'v  =  u'L'Lv  =  [0  ...  0,  /u'u  ]w  =  /l-r  •  w 
—  -  —  P 

The  joint  distribution  of  u  and  w  is,  therefore, 


(2.9) 


T(f+1) 

<2tt)  pr  (f-p+i) 


rs  -  (1 


-r)up}' 


f-p-i 

2 


dudw 


From  u  ,  transform  to  r  *  1  -  u'u  and  p-1  other  variables 

V  *2'  *•*'  *p-l  by 


(2.10) 


u.  =  (l-r)  cos$,  cos<J>_  ...  cos4>  ,  , 

1  12  p-1 

u.  =  (l-r)  1//2cos<J>_  cosd>_  ...  cos<j>  ,  sinif>  ... 

3  1  2  p-D  P-D+l 


( 3  —  2 ,  3  r  » *  *  p) 

Similarly,  transform  from  w  to  s  =  1  -  w'w  and  p-1  other  variables 

V  02 . 9p-l  ^ 


(2.11) 


V,  =  (1-s)  COS0,  COS0.  ...  COS0  ,  , 

1  12  p-1 

1/2 

v.  =  (1-s)  cos0,  cos0_  ...  cos0  ,  sin6  . , 

3  12  p-j  P-D+1 


( j=2,  3,  . . . , 
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The  Jacobian  of  transformation  from  u  to  r,  $  ,  is 

—  1  p-1 

.  j  p-1  p-2  p-i-1 

T-(l-r)  II  cos$. 

i-1  1 

and  a  similar  expression  in  s  and  6^  for  the  Jacobian  of  transformation 
from  w  to  s  and  the  0's.  Now  0  ,  and  d>  ,  vary  from  0  to  2n ,  the  other 

—  p-1  p-l 

0's  and  <fi 1  s  vary  from  -v/2  to  tt/2  while  r  and  s  vary  from  0  to  1.  Inte¬ 
grating  out  all  the  4>'s  and  all  8's  except  0  ,  we  obtain  the  joint 
distribution  of  r,  s  and  0  as 

fS£zl 

l2-12)  4~p-~lfFTfTp-;iT  tra  -  2  c°sP'26  drdsde 

where  0^  is  replaced  by  0 

The  joint  distribution  of  r,s  alone  can  now  be  obtained  by  integrating 
out  0  but  this  does  not  seem  to  yield  a  manageable  expression,  as  the 
bracket  in  (2.12)  will  have  to  be  expanded  in  a  series. 

3.  Moments  of  r,s  . 

Only  the  product  moment  of  r  and  s  is  difficult  to  obtain.  The 
mean  and  variance  of  r  (or  s)  can  be  very  easily  obtained  from  the 
marginal  distribution  of  r,  which  is  related  to  the  well-known  Hotelling's 
,  1 

T‘  by  r  =  — "■■■  ■  .  In  the  joint  distribution  of  u  and  v_  ,  given  by 

-W 

(2.3),  if  we  transform  to  z  *  (z^  ,  ...,  z  ]'  from  by 

(3.1)  v_  «■  (I  -  uu 1 )  1//2£  , 

we  shall  find  that  and  z  are  independently  distributed  as 

llE 

(3.2)  K (u |  f ) du  =  — — - lilZil  |i-uu'j2  du 

ttP/  (f-p)  j(f-p) 

(3.3)  and  K(£|f-l)d£  ,  respectively. 
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From  (2.2),  one  con  easily  show  that 
(3.4)  E(rh)  -  E(l-u'u)h  -  E|  X  -  uu'|h 

f(f+2h-p) 
(f-p) (f+2h) 


r(r?  *  hM?) 
r(f ♦  hK¥) 


This  will 
(3.5) 


also  be  the  h 
E(r) 


f-p+2 

f+2 


th 


moment  of  s  by  symmetry. 
v(r)  .  WV 


(f+2) 2 (f+4) 


This  leads  to 


Now  Cov(r,s)  ■  EUl-u’u)  (l-v'v)  }  -  E(r)E(s) 

-  E{(l-u'u) (l-z(I-uu')z)}  -  (E(r)}2  by  (3.1) 

-  E(r)  -  E{rl£'£-{£'u)  2] )  -  (E (r )  > 2 

(3.6)  =  E(r)  -  E(r)E(z'z)  +  E{r(z'u)2}  -  (E(r)}2  , 

as  z  and  r  are  independent.  Since  z  has  the  same  distribution  as  u 

with  f  changed  f-1  , 

E(z'z)  =  1  -  E ( 1  -  u'u)  with  f  replaced  by  f-1 


(3.7) 


=  -E_ 
f+1 


Hence  (3.6)  reduces  to 

(3.8)  Cov(r ,s)  =  ■  .~P^f~P+2^  -  +  E{r(z'u)2)  . 

(f+1) (f+2) 2 

Now 

(3.9)  E{r (£' u)  2 }  =  /  (1-u'u)  U'u)  2K(uj  f)K(£|  f-l)dud£ 

where  the  integration  is  over  the  range  of  values  of  u  and  £  such  that 
u'u  a.  1,  £'£  s  1  .  Transform  from  £  to  £  =  [£^,  . . .  ,  £p  1  by  the 
transformation 


£  =  L£  , 

where  L  is  already  defined  to  be  a  p*p  orthogonal  matrix,  whose  last  row 
is  u'/^u'u  .  Then, 

- -  1/2 

z'u  =  z'L'Lu  =  £'Lu  =  K  ►u'u  =  (1-r)  £ 

- -  -P -  P 
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Hence  (3.9)  reduces  to 

(3.10)  /r (l-r)K(u!  f )du  •  MjJ  K(ilf-l)d£*  K (r-r?)  •  •  *(£.'£>  ,  duo 


to  symmetry  of  the  distribution  of  ^  .  Now  £  hot  tho  inn  distribution 
as  u  with  f  replaced  by  f-1  and  hone*  finally*  (3.10)  reduce*  to 

pit-p+V _ L 

(f+4)  (f+2)  f*l 

The  covarianca  between  r  and  a*  therefore,  i»  (from  (3.5)) 

n  m  -2p(f-p*2) 

'  *  1  (f+l)  (f*2)<(f*4) 


Remarks » 

Wilks  considers  a  sample  of  size  n  and  has  a  Wishart  matrix  based 
on  n-1  degrees  of  freedom  as  deviations  are  from  the  sample  means.  He 
then  removes  two  observations  as  outliers  and  thus  his  (n-1)  -2  corresponds 
to  our  f  .  His  E(r)  agrees  with  our  result,  with  this  correspondence 
but  the  other  moments  are  in  error. 

Reference 

Wilks,  S.  S.  (1963).  "Multivariate  Statistical  Outliers,"  Sankhya,  Vol.  25, 
p.  407-426. 
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If  B  is  a  wishart  matrix  and  x,  £  are  two  vectors  of  p  components  each  having 
a  multinormal  distribution  and  if  all  these  quantities  are  independently  distributed* 
the  joint  distribution  of  the  two  statistics 


|b  ♦  h/  | 

| b  ♦  XX'  ♦  yy ' 


B  ♦  XX  ’ 


| B  ♦  XX '  ♦  | 


is  derived  in  this  paper.  The  correlation  between  r  and  s  is  also  obtained,  r  and  s 
are  related  to  Hotelling's  T?  and  are  useful  in  problems  of  testing  multivariate 
jutliers. 
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